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Abstract 

We study the phase space of the spherically symmetric solutions of Einstein-Maxwell-Gauss-Bonnet system 
nonminimally coupled to a scalar field and prove the existence of solutions with unusual asymptotics in 
addition to asymptotically flat ones. We also find new dyonic solutions of dilatonic Einstein-Maxwell theory. 
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1. Introduction 


In a recent paper [1] we have studied the global properties of the spherically symmetric solutions of 
Einstein-Gauss-Bonnet (GB) gravity nonminimally coupled to a scalar field. We refer to that paper for a 
full discussion of the relevance of GB gravity for physics and for a list of references. We only recall its role 
in the low-energy field theory expansion of string theory [2]. 

Our investigation showed that in four dimensions the only regular black hole solutions of the system are 
asymptotically flat. It is interesting to check if this property holds also for charged black holes. In fact, this 
is not the case. Already in absence of GB coupling some degenerate charged solutions of dilatonic gravity 
with unusual asymptotic behaviour are known [3]. Actually, analytical solution to dilatonic gravity with 
nontrivial Maxwell field are well known in absence of GB coupling [4,5]. If only magnetic or electric charges 
are present, the solutions display only one horizon and are analogous to Schwarzschild black holes. If instead 
both charges are nonvanishing, they possess an inner and an outer horizon and the causal structure is similar 
to that of the Reissner-Nordstrom solution. Asymptotically, they are flat except in the degenerate cases 
cited above [3]. 

In the case of dilatonic Einstein-GB models, magnetically charged solutions cannot be obtained in 
analytical form, but have been studied perturbatively [6] and numerically [7,8]. Only asymptotically flat 
solutions have been investigated in this context. The existence of solutions with different asymptotics 
can however be decided through the study of the global properties of the phase space associated with the 
dynamical system, as has been shown in the literature in analogous situations [9]. It must be noticed that 
the search for solution is made easier in the Einstein case by the existence of a duality between electrically 
and magnetically charged solutions, that is no longer available in the presence of GB corrections. 

Also interesting is the study of the emergence of branch singularities in the metrics. This is a common 
feature of GB models [10,7], and can give rise to naked singularities. This topic can be investigated by means 
of an expansion of the fields near the outer horizon [11]. 

In this paper, we study the dilatonic Einstein-GB system coupled to Maxwell field, and classify the 
solutions by their asymptotic behaviour . We find that their behaviour is the same as in the dilatonic Einstein 
case, except for some non-asymptotically flat purely magnetic solutions, which do not have analogues if GB 
terms are included. In the course of the investigation, we also find new dyonic solutions of dilatonic Einstein 
gravity with unusual asymptotics. We also study the near-horizon expansion of the metric functions, finding 
the conditions for the absence of naked branch singularities and show the existence of solutions with multiple 
horizons. 

2. The dynamical system 

In this section, using the methods of ref. [1], we derive the dynamical system associated to spherically 
symmetric solutions of the field equations of dilatonic Einstein-GB-Maxwell theory. For further technical 
details we refer to that paper. 
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The four-dimensional action of four-dimensional Einstein-GB gravity coupled to electromagnetism is 


I = J d*x [7^ - 2(a^)2 -K ae - F^)] , (1) 

where S — + TZ"^ is the GB term, F is the Maxwell field, ^ is a scalar field and a is 

a constant with dimension (length)^. This action also emerges in the low-energy expansion of string theory 
[2], and in this case the constant a is proportional to the slope parameter. Note, however, that in string 
theory also an axion field may couple to the gravitational and Maxwell fields, and this may change notably 
the properties of the solutions [ 12 ]. 

For the metric we adopt the spherically symmetric ansatz 

ds"^ =—Adt"^ + A~^ dr^ + , ( 2 ) 

where A and R are functions of r. In these coordinates the generalized Maxwell equations 

= 0 (3) 

admit the spherically symmetric solution 

g20 

F = Qe -tpt dt a dr + sin 9 dO A dcj), (4) 

R^ 

where qe and qm are the electric and the magnetic charge, respectively. 

The field equations for the gravitational and scalar degrees of freedom following from (1) and (4) are 
more conveniently written in terms of new variables such that the metric and the scalar take the form [ 1 ] 

ds^ = -e^+^dt‘^ + e^‘^-'^-^d^‘^ + e^^-^-^dn^, 2(j> = tj - x, (5) 

where rj, x and A are functions of They read 

C" - - x'^)iv' + x' - 2C')e^^"^‘^]' = -f a{r]'^ - - {rj' + x' - 2C')^e^^“'‘^], 

rj" + a[Ar]'- {rj' + x' - - x'^ + v'x' - = 2aq'^e'^^, 

x" - a[4x'e^^"^‘^ -k {ri' + x' - 2C')(??'^ - 2x'^ - ??'x' + 2r]'= 

2aql,e^^ + airC^ - x'^)[4e^^“^^ - (ry' -k x' - 

2^ - - x'^)[4e^^”^‘^ - 8 ( 77 ' + x - 2Cfe^^-^^] = 2e^^ - 2aql,e^^ - 2aqle^'^ . ( 6 ) 

It is interesting to notice that in spite of the presence of the higher-derivative GB term, the field equations 
are second order. 

We can then write ( 6 ) in the form of a first order dynamical system by defining new variables 

As: = C', E = ry', IE = x', Z = , U = V^e^, V = V^e^, (7) 
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which satisfy the differential equations 


Z' = XZ, U' = YU, V' = WV. (8) 

In terms of these variables, the field equations become 

f ' r 1 

X' - <^ - W'^){Y + W- 2X)^ I =Z^ + (Y^- W^) -(Y + W- 2Xfj (9) 

r 1 V 

AYZ"^ -{Y + W- 2X){2Y^ - + YW - 2XY)^ ^ = 2qlu'^, (10) 

f r 1 t/2 

W' - j l^iWZ^ + (Y + W- 2X){Y^ - 2W‘^ - YW + 2XW}^ 

+ {y^ - W^) -{Y + W- 2Xf] (11) 

L \ Ai 

subject to the constraint 

r 1 

E = 2X^-Y"^ -W"^ - 2Z‘^ + 2qly‘^ + 2qlu'^ - (F^ _ vh^) ^4Z^ _ 3(y + ty - 2Xfj (^2) 

In general, the dynamical system is six-dimensional, but for purely magnetic solutions {qe = 0), the variable 
U is ignorable and does not appear in the system, while (10) yields the first integral 

r 1 y^ 

Y + \4YZ^ -{Y + W- 2X){2Y^ - W^ + YW - 2XY)^ ^ = «■ (13) 

In this case the dynamical system is effectively four-dimensional. We shall nevertheless work also in this case 
with the full six-dimensional phase space, because this facilitates the discussion of the critical points. For the 
same reason, we shall not use the constraint (12) to eliminate V'^/Z'^ and reduce by one the dimensionality 
of the phase space. 

The system can be put in canonical form solving (9)-(ll) for the first derivatives, as 

X' = F{X,Y,W,Z,U,V), Y' = G{X,Y,W,Z,U,V), W' = H{X,Y,W, Z,U,V). (14) 

The expressions so obtained are awkward and shall not be reported here. 

3. Charged solutions of dilatonic Einstein gravity 

In the following, the knowledge of the exact solutions in absence of the GB term will prove useful. These 
were first discussed in [4], and later some degenerate cases with nontrivial asymptotics were found in [2] for 
the magnetic monopole. We shall generalize these solutions to the dyonic case. An important property of 
the dilatonic Einstein-Maxwell theory is the presence of a duality invariance for the interchange of q^ and 
qm, with (j) —(j). 

In absence of Gauss-Bonnet contributions, the field equations (6) reduce to 

v" = 2qy\ X" = 2qle^\ (15) 
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subject to the constraint 


2('2 _ ^/2 _ ^/2 _ 2 e 2 C ^ 2 ( 7 ^+ 2 qy ^ = 0 , 


(16) 


where we have absorbed a into the definition of the charges. In these variables, the duality invariance 
corresponds to the interchange of r] and y. The integration of (15) yields 


= 


2ae“f 


V2 


72 


ce 




(17) 


1 - e2< ’ \qe\ 1 - e2K?-«i) ’ \q^\ 1 - e 2 c(«-« 2 ) ’ 

where a, b, c, and ^2 are integration constants and the origin of ^ has been fixed arbitrarily. Substituting 
in the constraint (16), one gets 2a^ — = 0. If one also requires that the solution has a regular horizon, 

one must impose that is regular at the horizon, ^ ^ — 00 , which implies 2a — h — c = 0. 

Combined with the previous constraint, this yields a = b = c. Substituting these results in (5) and defining 
a new coordinate r = J e'^‘= = 2a/(I — e2“?), one obtains the solution 


r(r - Th) 


ds^ = — 7 N / N 

[r + re)[r + rm) 
- 24 , ^ 7 + re) 

{r + Tm) ’ 


^, 2 +(’- + r.)(r + r„) 


r(r - Th) 


.){r 


^) dn^ 


(18) 


where rh = 2a and Vg and are functions of a, and ^2 • This solution was first discussed in ref. [4]. The 
mass M and the charges of the black hole are related to the parameters of the solution by 


2M = Th + Te + Tr, 


2 _ re{rh + re) 


2 r^{T}^ “f r^) 

dm 


(19) 


2 ’ 2 

The metric (18) has two horizons aX r = Vh and r = 0, a singularity at r = max(—rg, —r^), and is 
asymptotically flat. The special case when ge = 0 was also derived in [5] and corresponds to re ^ 0. 

Some degenerate solutions with unusual asymptotics arise when or ^2 vanish. For example, for ^2 = 0, 
one has 


2 r{r-rh) 


ds = — - 


J^2 , T+Te j 2 , r , ^ jo2 -24 ^ + 

dt H—;- -dr +{r + re)dil <= ^ — 


( 20 ) 


r + re r(r - th) ’ 2g^ 

For = 0, the metric has the same form, with Vg r^, while the scalar becomes = 2q/./{r + r^,). In 
different coordinates the metric (20) can be also written 


ds^ = — 


{R^ - rg){R^ - ro) 


dr 






( 21 ) 


i?2 (i?2 - re)(i?2 _ ro) 

with = i?^/2g^, and ro = rt + Vg. The solution (20) generalizes that found in [3] for the purely 
magnetic case, which is recovered in the limit rg —> 0. In the general case, the metric (20) displays an 
outer horizon at r = r?i, an inner horizon at r = 0, and a singularity at r = —rg. Its asymptotic behaviour 
is intermediate between Schwarzschild and anti-de Sitter and its causal structure is similar to that of the 
Reissner-Nordstrom metric. 

For = ^2 = 0, one obtains a previously unknown solution, of the form 




ds^ = -- 


‘^dedm 


‘2qgqm 
r(r - Th) 


dr'^ + 2qeqmdn‘^, 


Qm 


( 22 ) 
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This metric is the direct product of two-dimensional anti-de Sitter spacetime with a two-dimensional sphere 
and is therefore of the Bertotti-Robinson form. It has a horizons at r = 0 and r = Vh, but no singularity, 
and both nonvanishing electric and magnetic charge. 

In all cases, extremal black holes are obtained in the limit a = 0, corresponding to r/, = 0, or, in terms 
of the physical parameters, 2M^ = (^e ± 

Finally, we notice that in terms of the variables of the dynamical system the solution (17) reads 

X = acoth(a^), F = acoth[a(^ — ^i)], IF = acoth[a(^ — ^ 2 )], 

^ _ 

sinh(aC)’ \/2|9e| sinh[a(^ - ^ 1 )] ’ sinh[a(^ - ^ 2 )] ’ 

with the horizon at ^ ^ — 00 , and spatial infinity at ^ = 0. From the explicit knowledge of the solutions or 
by standard methods, it is then easy to construct the phase space portrait. Phase space is invariant under 
Z —Z, U —U or F ^ —V, thus we shall limit our study to positive values of these variables. 

The physical trajectories lie on the five-dimensional hyperplane E = 0. The critical point at finite 
distance lie on the surface Zq = Uq = Vq = Pq = 0, where 

= 2^2 - f2 _ (^23) 


but only points with Xq = Yq = Wq correspond to regular horizons. The eigenvalues of the linearized 
equations around the critical points are 0(3), Xq, Yq, IFq. 

The phase space at infinity can be studied defining new variables 


Y 

IF 

z 

u 

V 

— 

w = —, 

Z '=■ — 

u = —, 

V =■ — 

X’ 

X ’ 

X’ 

X’ 



(24) 


The critical points at infinity are found at to = 0 and 

a) Mo = Wo = -Zo = 0, y = yo, w = wq, with Vq+Wq = 2. These are endpoints of trajectories lying on 
the surface at infinity. 

b) uq = vq = 0, zo = 1, yo = yJo = 0. These are the endpoints of the asymptotically fiat solutions (18). 

c) Mo = 0, Wo = XjZo = I, yo = 0, wo = 1. These are the endpoints of the solutions (20). 

d) Wo = 0, Mo = 1 /a/^, Zo = 1, 2/0 = 1) Wo = 0. These are the endpoints of the dual solutions to (20). 

e) Mo = wo = 00 = 1, 2/0 = iwo = 1. These are the endpoints of the solutions (22). 

The eigenvalues of the linearized equations are 

a) 0(3), 1, yo, iwo- 

b) -1(5), -2. 

c) 1, -1(3), -2(2). 

d) 1, -1(3), -2(2). 

e) 1(2), -1, -2(3). 


6 



In absence of electric charge, the system reduces to a four-dimensional system with U = Y = 0, and similarly 
for qm = 0, it reduces to a system with V = W = 0. 

The phase space portrait is the following: regular solutions start at Zq = Uq = Vq = 0, Xq = Yq = 
Wo = a, for some value of the parameter a. At infinity, a three-dimensional family of trajectories is attracted 
by points b) and a two-dimensional bunch by points c) and d). Finally, a one-dimensional bunch ends at the 
critical point e). 

4. Phase space of the Einstein-GB system 

We discuss now the phase space associated with the Einstein-GB system (8), (14). Since E in (12) 
is a first integral of the dynamical system, all the trajectories are confined to surfaces of constant E, but 
only those with vanishing E are physical. We shall therefore consider only the critical points lying on the 
hypersurface E = 0. The discussion of the phase space is complicated by the presence of a singularity in the 
field equations at Z = 0. Since some of the critical points of the dynamical system lie on this hypersurface, 
one must be careful in taking the limits when approaching these points. In particular, when Z 0, also 
E —> 0. Although the duality invariance is no longer valid, the phase space is still invariant under Z —>■ —Z, 
U —U or V ^ —V, as in the pure Einstein case. We can therefore consider only positive values of Z, U 
and V. 

The critical points at hnite distance are defined as the points were the r.h.s. of eqs. (8) and (14) vanishes, 
and in analogy with the Einstein case are found to lie on the surface Zq = Uo = Vo = Po = 0, with P given 
by (23). We are interested in regular black hole solutions. A regular horizon is present if is 

finite at the critical points. Moreover, the dilaton is regular at the horizon if e~‘^^ = is finite there. 

These conditions are equivalent to the requirement that at the horizon ~ ~ Since the critical 

points correspond to ^ ^ ±oo, with Z^ ~ ~ ~ ^ one obtains that the horizon is 

regular iff Xo = lb = Bb = a, for some real parameter a. To simplify the following discussion, we shall only 
consider points that satisfy this constraint. 

The nature of the critical points can be investigated linearizing (8) and (14) around them. One obtains 
6Y' = SX' = 6W' = 0, SZ' = a 6Z, SU' = a SU, SV = a 5V, 

with eigenvalues 0 (3), a (3). It results therefore that the critical point is an attractor for a < 0 and a 
repellor for a > 0. 

The asymptotic behaviour of the solutions is dictated by the nature of the critical points at inhnity 
[8,1]. These can be investigated using the variables (24). The field equations at infinity {t 0) read 

i = - ft, y = 9 - fy, w = h- fw, 

i = (1 - f)z, u={y- f )u V = {w- f)v, 

where a dot denotes td/d^, f = t‘^P{l,y,w, z,u,v), g = t^G{\,y,w,z,u,v), h = t^P[{\,y,w, z,u,v), and the 
variables are subject to the constraint 

2 - y"^ - - 2z‘^ + 2ql^v^ + 2qlu^ - (j/ - w^) Uz^ - 3(y + w - 2)^1 ^ = 0. 
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The critical points at infinity are placed at to = 0, with 
a,) uo = Vo = zo = 0, y = yo, w = wq, with y'^ + = 2. 

b) uo = vo = 0, Zo = 1, yo = Wo = 0. 

c) Vo = 0,uo = 1/ Zo = l,yo = 1, wo = 0. 

^0 = ^^yo = wo = 1 . 

The points a), b) and c) coincide with those obtained in the Einstein case, while the point e) is displaced due 
to the GB contribution. The most interesting feature is however that the critical point d) has disappeared. 
In particular, no solution with vanishing electric charge and unusual asymptotics exists in the GB case. 

The eigenvalues of the linearized equations are 

a) 0(3), 1, yo, wq- 

b) -1(5), -2. 

c) 1, -1(3), -2(2). 

e) 1, -1, -2(2), -i(l ± ^9+16/ql). 

They have the same sign as in the Einstein case, and hence analogous properties. 

The resulting phase space portrait is rather similar to that of the Einstein system, except for the critical 
points d), that in the GB case are absent. In particular, all the trajectories corresponding to solutions with 
regular horizons start from the critical point at Zq = C/q = Vq = 0, Xq = Iq = ITo = a > 0 for a given value 
of the parameter a, and end at the points b), c) or e). 

From the location of the critical points at infinity, one can deduce the asymptotic behaviour of the 
solutions [9,1]. The behaviour coincides in all cases with that found in the Einstein limit. In particular, 
in case b) solutions are asymptotically flat with asymptotically constant scalar field, while in case c) they 
behave like (20) and in case e) like (22). The existence of solutions with a = 0 also suggests the existence 
of metrics with multiple horizons, i.e. with multiple roots of the metric functions, that should correspond to 
extremal black holes. 

5. Branch singularities 

In the previous section, we have studied the existence of solutions with regular horizons. Our method, 
however, does not give a complete characterization of the black holes. In particular, it does not permit to 
study the region inside the horizon and to ascertain the existence of inner horizons^ Moreover, we are not 
able to relate the parameters a, Qm and qe to the radius of the horizon. On the other hand, in GB models, 
branch singularities usually appear [10,7]: the requirement that they are shielded by an horizon imposes a 
minimal value for the horizon radius (and consequently for the mass of the black hole). In order to investigate 
the presence of these singularities, it is useful to study the behaviour of the fields near the outer horizon 
[11]. The calculation is most easily performed in coordinates of the form 

ds^ = + R^d^f. (25) 


^ Some results on this topic have been discussed using numerical methods in [7]. 



In this gauge, the field equations read 

, I'm' - 2 \ 

^ +^+sr = 


2«e \ _ 


i?2 

A' 1 - 

'r W~ 


- + (1 - e-^) + |(a*' - A')jJ - a{qi,e-^^ - (26) 

= [(1 - - 2(1 - e-^){r - 20'^)] + a{qle-^^ + (27) 


^ ^ ^‘^(l-3e ^)fi'<j)'- aiql^e ^'*’+ qje^'*’)^, 


8ae-2<^ 


e-" U" + ^(m'-3A'))<^' + m'(<^"- 2<(.'^) +2a(g^e-2^ + g,V^)|-. (29) 


- 2())'2 _ 


In order to find an analytical expression for the minimal value of the radius, we assume that near the 
horizon R = Rh the fields behave as 


e-^ = A(i?-i7^) + ..., = M(i?-!?,,) + ..., 

(j) = 4‘h + ~ ^h) + . • . , 


for some constants A, M, and Let us first consider the case qe = 0. Substituting into the field 
equations , one obtains from (28) 

(30) 


Rh 1 - / ’ 


where 


Rl ’ 




Moreover, combining (26) and (29), it results that 


8a2e-4<^'* 


A + ( 1 + RRhf — 


Aae ^‘^'“< 7,1 


RRl ’ 


and hence, substituting (30) in (32), after some algebraic manipulations one obtains a second degree algebraic 
equation for /: 

[1 — (1 + a + a^)gm]/^ — [1 — (1 + 2d + a^)gm + + 2 [3A^ ~ 2d(l + 5a)gm + = 0, (33) 

where we have defined a = 

In order to have real solutions for /, the discriminant A of the equation must be nonnegative, i.e. 

A = (1 — gmY[l — 6d2 + 2d^(l + 3d + 3d^)gm + d'^g^] > 0. (34) 

The inequality holds for any g^ as long as d^ < 1/6, i.e. R^ > /A^/Q. The same restriction holds 

also for neutral solutions [11,7]. However, in the charged case (34) can be satisfied also when d^ > 1/6, if 


-(2 + 2d + 3d2) + a/ 3(1 + 4d + 9d2 + 6d3 + 3d4) 
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The same result was also obtained in [8] for SU(2) Yang-Mills monopoles. It is also interesting to remark 
that the vanishing of the prefactor in (34), = 1, corresponds to extremal black holes. 

When Qe ^ 0, one obtains similar results. In this case, defining ge = q^/Rl, 

A = (I - 5m - 5e)^[l - 6 d^ -b 2 d^(l -b 5d'^){gm, + ge) +d‘^i9m+ 9e)'^ + 6 ( 5 ^ “ 9e)], (35) 

which is nonnegative for d^ < 1/6. If instead d^ > 1/6 a complicated constraint must be obeyed by the 
charges. Again, the case gm + 5e = 1 corresponds to the extremal limit. 

6. Conclusions 

The study of charged solutions of the dilatonic Einstein-GB model has revealed the presence of some 
difference from the pure Einstein case. In particular, non-asymptotically flat solutions with vanishing electric 
charge are absent in the GB case. The conditions for the absence of naked branch singularity is identical 
to that of the neutral case, but the possibility exists of solutions violating this condition, if the charges 
satisfy some inequalities. Also interesting is the existence of extremal solutions, which were not considered 
in previous investigations. 
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